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We analyze the spectrum of the electric-current autocorrelation function (noise power) in the 
Anderson impurity model biased by a finite transport voltage. Special emphasis is placed on the 
interplay of non-equilibrium effects and electron-electron interactions. Analytic results are presented 
for a perturbation expansion in the interaction strength U. Compared to the non-interacting setup 
we find a suppression of noise for finite frequencies in equilibrium and an amplification in non- 
equilibrium. Furthermore, we use a diagrammatic resummation scheme to obtain non-perturbative 
results in the regime of intermediate U. At finite voltage, the noise spectrum shows sharp peaks at 
positions related to the Kondo temperature instead of the voltage. 

PACS numbers: 85.65. +h, 73.63.-b, 63.22.-m 



I. INTRODUCTION 

The Anderson impurity model (AIM) is probably the 
best studied interacting quantum impurity systemP^l 
Nonetheless, there are still a number of open ques- 
tions, especially concerning the physical properties of 
its non-equilibrium multi-terminal incarnations. These 
are archetypical models to describe transport proper- 
ties of nanoscale quantum dots. In recent years not 
only the classical transport characteristics, such as lin- 
ear conductance and full current-voltage relations, have 
been studied, but also more complex quantities, like the 
full counting statistics, have been explored. Thereby a 
large variety of different analytical as well as numeri- 
cal techniques have been applied to the AIM (see for 
instance [5lU8|). However, most of the studies have con- 
centrated on static or integral energy independent prop- 
erties. Aiming at the ultimate practical applications of 
the devices in question, the knowledge of their frequency- 
dependent noise would undoubtedly be of great advan- 
tage. Currently, research in that direction appears to 
gain momentumliSHU] The frequency-dependent noise is 
not of mere theoretical interest as even the frequency- 
dependent current cumulants of higher order are by now 
measurable in experiments.^ 

Depending on its parameters the AIM shows many in- 
teresting phenomena, mostly with signatures in transport 
properties of the two-terminal setup. Probably the most 
prominent one is the Kondo effect which is generally only 
visible for large on-site repulsion U. Very few analytical 
methods starting from the proper AIM are able to ac- 
count for this correlation effect, as the analytical solu- 
tion for the AIM at arbitrary U in non-equilibrium does 
not exist. As far as the integral transport characteristics 
are concerned, it was shown in Refs. |9|ll)23j that the 
current-voltage relation as well as the full counting statis- 
tics can be calculated up to the cubic order in applied 
voltage and for arbitrary U . For the frequency-dependent 
noise the non-interacting situation [resonant level model 
(RLM)] in non-equilibrium was addressed in Refs. [19124) . 



A number of interesting frequency-dependent noise re- 
sults were obtained by use of the equation-of-motion 
techniquepS and Filippone et al. recently investigated 
similar quantities in a slightly different context!^ Nev- 
ertheless, the general crossover from weak to strong U 
— eventually entering the Kondo regime — was not yet 
addressed for the frequency-dependent noise. With this 
work we want to close this gap and provide results valid 
for arbitrary bias voltage and hybridization parameter. 

The paper is organized as follows. In Sec. [Tljwe start 
with a simple perturbative expansion in U. Using the 
recently derived analytical results for the non-equilibrium 
self-energy we analyze the interplay of finite voltage and 
on-site repulsion and identify the related features in the 
noise spectrum. In Sec. |III| we derive a diagrammatic 
resummation scheme (which is similar to the one used in 
Ref. [37] for the shot noise) that enables us to generate 
non-perturbative results. We find that for larger U the 
noise spectra show distinct features which can be related 
to the Kondo effect. Section |TV| concludes the paper with 
a summary of the results. 



II. PERTURBATIVE EXPANSION IN U 

We model the system by the two-terminal AIM Hamil- 
tonian 

H = H L + H R + H d + H T + H! (1) 
H d = J2 A( & d ° 

a 

H T = 5>(4^+4^r+H.C.) 

u 

Hz = [/4(i)d t (i)4(iW)- 

Hd describes the spin(cr) degenerate dot level at energy A 
with the fcrmionic creation and annihilation operators d\ 
and d a . Ht is the tunneling Hamiltonian with tunneling 
amplitude 7 to the left (L) and right (R) leads^ and 
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Hi models Coulomb interaction on the dot. Hl (Hr) 
describe free electrons in the left (right) lead, held at 
constant chemical potential [i L = V/2 (fi R = —V/2). 
The principal quantity to calculate is the current noise 
spectrum 



S(u>) = f dte-^((I(t)I(0))-(iy 



(2) 



given as the Fourier transform of the current autocor- 
relation function. As was pointed out in Refs. [2"5H3"U] 
this unsymmetrized definition of S(uj) is the most nat- 
ural one for quantum systems: negative frequencies de- 
scribe the spectrum of fluctuations which emit energy 
to the environment, positive frequencies are related to 
the absorption spectrum. Furthermore, in order to be as 
close to the experimental situation as possible, we define 
the total current operator through the constriction as the 
average through the two contacts (we use units in which 
e = H = c = 1 throughout), 

I ( t )= i l E [4(*K(*)--Rt(*K(*)-H.c] . (3) 

As a first step, we evaluate the current correlations 
(I(t)I(t')) in the non-interacting case (U — 0), when 
the spin sector completely factorizes. The noise can be 
rewritten in terms of Keldysh Green's functions (GF) 
Gf tap {t- f) = -i(T c a(t)ft{t')} where a,fi are lead 
{L,R) or dot (d) creation or annihilation operators, k 
and I are Keldysh indices of the times t and t' and Tq is 
the Keldysh time ordering operator. The average ( ■ )o is 
over the non-interacting part of the system only for which 
the GF are known exactly. In particular, w e defin e the 
dot GF by Df = G%\ 



^ k , l dd for which one obtains^ 



D a (u) 



W-A+2irln(u>)-l| 
(w-A) 2 +4r 2 
2ir[n(w)-2] 

(w-A) 2 +4r 2 



2iTn(ui) 

(w-A) 2 +4r 2 

-w+A+2ir[n(aj)-l] 

(w-A) 2 +4r 2 



(4) 



Here we defined n(cj) = 9(V/2 - oj) + 6{-V/2 - w) with 
the unit-step function 6. The electrodes are described in 
the framework of the wide flat-band model with constant 
density of states po- Furthermore, we use 2T = 2wy' 2 po 
as the unit of energy if T does not occur explicitly. In the 
following we restrict the analysis to the zero-temperature 
case as one expects only a smear-out effect in the noise 
spectrum at finite temperature. Within these approxima- 
tions we obtain the following analytical result for U = 0, 



77 <7,a' = ± I 



arctan (V/2 + a' A + ulj) + arctan (V/2 + a' A) 
(V/2 + a' A + auj) 2 + 1 



■P{u,oV) In 



(V/2 + a' A) 2 + 1 



(5) 



Here P(x, y) = 9(x) — 9(y + x) is a 'window' function and 
Q(x, y) = 0(x) + 9(y + x) is a double-step function. Our 
result (|5J) is in agreement with the numerical results of 
Ref. |19[ . It can be derived from the general expression 
for a two-terminal system with the Lorentzian-shaped en- 
ergy dependent transmission coefficient (see e. g. [24] ) 

T(u) =4r 2 /((w- A) 2 +4r 2 ). 

The generic features of the non-interacting noise spectra 
are characteristic steps in the unsymmetrized noise spec- 
trum (and cusps in the symmetrized version) at energies 
corresponding to ±A ± V/2. 
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FIG. 1: Unsymmetrized noise spectrum for the non- 
interacting case. The negative/positive ui part can be in- 
terpreted as emission/absorption of energy by the system 
to/from the environment. Steps emerge whenever the energy 
uj of a fluctuation is about ±V/2 ± A. The inset shows the 
symmetrized spectrum for the same parameters. 

It is well known that the most prominent feature of 
the Kondo effect is the emergence of a sharp resonance 
around the Fermi edge of the metallic continua coupled to 
an impurity. An effective model at energies much smaller 
than the Kondo scale Tk is a resonant level mod el, where 
the role of hybridization V is given by T r * 2 * 32 1 There- 
fore our result ^ can formally be used to model the 
behaviour of the noise in the deep Kondo limit. Indeed, 
an explicit fit of the data from [35] shows a very good 
agreement with the experiment provided that effective 
T w 4.9/ieV. 

Now we turn to the perturbative expansion in U. As 
long as one is interested in the electron-hole symmet- 
ric case A = —U/2 there is no linear contribution in 
U and the expansion starts with terms of quadratic 
orderP^l We restrict the following analysis to this spe- 
cial case since it contains most of the interesting fea- 
tures. In contrast to the interaction-free case there is 
no spin-sector factorization anymore. The expression 
for the noise can still be obtained in a straightforward 
way. It contains a large number of terms which can 
be subdivided into two large classes. We discuss exem- 
plarily at full extent one representative term which con- 
tains dot operators and those for the left electrode. All 
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other terms can be treated equivalently. The two classes 
of terms are the parallel-spin contribution S dLdL — 
{d\{t)L^(t)d\(0)L^{0)) and the anti-parallel -spin contri- 
bution S>jL dL = (dl(t)L-t(t)d^(0)L±(0)) . In comparison, 
both have very different structures and thus need to be 
considered separately. 





(3) 




O 



(4) 



FIG. 2: Diagrammatic representation of the parallel-spin cur- 
rent correlations in the second order perturbation theory for 
the electron-hole symmetric case. All arrows represent non- 
equilibrium Keldysh GF exact in the U = model. The 
wiggly line represents instantaneous Coulomb interaction. 
Dashed green arrows describe particles with opposite spin 
than solid black ones, t and label the two external points 
of the currents in equation [2] Information about Keldysh 
indices and lead affiliation (especially if of homogeneous or 
inhomogeneous kind) of the GF are omitted. 

For the parallel-spin contribution in second order we 
obtain 



£ « / 



7 f duii 
dwi du> 2 du> 3 



~4 I ^ G 0,Ldi Ul ) G 0,Ld} 



k,l=± 



(2tt) s 



{B 1 + B 2 + B 3 - B 4 ) 



(6) 



with the definitions 



B l =D l q ^u> x )D«\w 1 + u 2 -uj 3 )D 1 «{u 2 )D«-{uj 3 ) 
x G old( w 3 - w)G+' (u; 3 ) , 



B 2 =D l k (u) 1 )D^(ij 1 +uj 2 - Lj 3 )D l k (u> 2 )D k+ (ij 3 ) 

x G (tZd( w 3+w)G , o ! L(w 3 ), 
B 3 =D lk {ui)D h \u 1 +uj 2 - lu 3 )D 1 -{lo 2 )D% + (lj 3 - u) 

X G (tL( W 3)Go^(^2-w), 

B± =D lk {u 2 )D k \uj 2 - uj)D 1 +{u 3 - w)Dj-(wi) 
*G+> Ld {u 3 )G- k Ld {ui-u). 

The corresponding diagrammatic representation is shown 
in Fig. [2] While the diagrams (1) and (2) can conve- 



niently be reduced to the self-energy of the second or- 
der, for which analytic expressions exist even in non- 
equilibriumJSl the evaluation of the other two diagrams 
involves an extra energy integration, which can conve- 
niently be performed numerically. 

No reduction in terms of the self-energy is possible for 
the anti-parallel-spin configuration contributions given 
by 



S. 



n 

dLdL 



M = -u 



2 7 



. 1 L J 



k,l=± 



dwi dcu 2 du 3 

(iff) 1 



{D 2 + D 3 ) 



(7) 



Here 



D 2 =D l ^(ui)D^(ui +lo 2 - w 3 )D l Q -{io 2 )D k+ {u 3 - u) 

x G+^MG^K - w) , 
D 3 =D l k {uj 1 )D k -{u Jl +uo 2 - uj 3 )D lk {u 2 )D k+ {u 3 - u) 

x G^ Ld ((Ji +uj 2 - uj 3 )G+ l Ld (uj 3 ) . 

which can be represented by the diagrams shown in 
Fig. [3] It turns out that all diagrams contribute to the 
noise power with the same order of magnitude. 
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FIG. 3: Diagrammatic representation of the anti-parallel-spin 
current correlations in the second order perturbation theory 
for the electron-hole symmetric case. 

An important benchmark for the full result is the lim- 
iting behavior at zero-frequency [the shot noise 5(0)] for 
the second order in U corre ction , for which analytical 
results exist for small voltages^H 



3t^4 



5(0) = 6U'\V\ 6 /tt 6 T 



(8) 



This limit is perfectly reproduced by our calculation of 
the full noise up to voltages V/Y sa 1. Note that it was 
shown in Ref. 19J that S{uj) is a purely real quantity 
which is also reflected in our calculation. 

Figure [4] shows the second order correction to the full 
noise as a function of frequency for different values of V. 
The overall scale of this contribution is about 1% of the 
U = value for U/Y = 2 and scales with U 2 . As U = 
is a fixed point in the RG-sense, it should be valid for 
small U. In equilibrium a general suppression of noise 
is observed (see solid curve in Fig. [4]). Since the dot 
can be considered to have a finite capacitance one is ef- 
fectively confronted with a high-pass filter, which leads 
to the above mentioned effect. The situation is com- 
pletely different in case of a finite bias. First of all, even 
a very small V leads to non- vanishing noise around u> = 0, 
where the noise is zero in equilibrium. This remarkable 
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FIG. 4: Correction to the noise spectrum in the second order 
perturbation theory in the resonant configuration A = —U/2 
in units of U 2 . 



behavior can be traced back to the imaginary part of the 
self-energy, which is algebraic in uj in equilibrium but 
immediately acquires a constant ~ V term after a bias 
is switched on. As this part of the self-energy describes 
dissipation, it is natural that V ^ leads to enhanced 
fluctuations in the low-energy sector. The form of this 
contribution is nearly Lorentzian with a width ~ T, as 
opposed to the simple step in the non-interacting case. 
For V > T this peak splits into two — similar to the 
doubling of the step in the non-interacting case — and 
as the voltage grows both peaks shift roughly linearly in 
V. From the mathematical point of view this behavior is 
understandable as the noise correction is composed out of 
GF with Lorentzian shape that are cut off by the chemi- 
cal potentials in the leads. This step (non-interacting) vs 
peak (interacting) behavior of the noise should be clearly 
distinguishable in the relevant experiments. 



III. DIAGRAMMATIC RESUMMATION 
SCHEME 

Although the perturbative expansion up to U 2 might 
yield reasonable results^ — at least for large voltages — 
this is not the generic situation encountered in experi- 
ments (see e. g. Ref. [311 )■ The same applies to the oppo- 
site limit of infinitely strong interactions for which a num- 
ber of interesting results became available recentlyPS. In 
the following, we address the range of intermediate U 
by using a diagram resummation technique, which was 
presented in a simpler version earlier in Ref. [27] . 

The full noise can be expressed with help of a two- 
particle GF (which we call T klmn , with 4 Keldysh in- 
dices), 



-^klran 



{tiMMM) = (T c dM)dl{t l 2 )d a ,{tf)dl,{tl)) 



S{oj) = 



FIG. 5: Reduction of the noise to the 2-particle GF Y. A 
distinction between correlations of anti-parallel-spin particles 
and of parallel-spin particles has to be made but was omitted 
in this figure. To express the noise in this way one needs two 
free lead GF Ql/r expressed by the dotted arrows. 



which involves only dot operators. Therefore, two free 
lead GFs are separated from the average values, see Fig. 
[5] Then the frequency-dependent noise can be expressed 
as 



S(w) 



z E E « 



dwi dw2 



k,l=± cr,<T'=t",4. 

r^r( Wl , w + Wl ,w 2 ) 5 l fc (- Wl ). 9 z'(- W2 ) 

r^+7 Z (wi,cj + u)i,u 2 )g± k (-uji)g l r (u) - lo 2 ) 

- T+ k Jr (u 1 + u,u 1 , u 2 )g- + ( - wi )gZ l (-oj 2 ) 



- E 

4 ^ 

fc=±=-i 



dwi 

~27T~ 



D^i-Ms+M- (10) 



kl 1 



Here we have taken the four dimensional Fourier trans- 
form of r WTO "(t 1 ,t 2 ,t3,t 4 ) and neglected its 4 th argu- 
ment, which is possible because of energy conservation. 
Furthermore, we have used the definition 

gf{u)= 7 2 [ 5 £'M±^H] (11) 

where the left/right lead free GF are denoted by g^/R- 



v x <^ x 

ivl-u'i }+,u>7+v' Q+ 



(9) 



FIG. 6: The complete summation over all diagrams that con- 
tribute to r in our approximation. All arrows in this case 
represent the dot GF Do of the U — model which is exact 
in tunneling. Every second diagram is for anti-parallel-spin 
current correlations. 
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In general, the perturbation expansion for this kind 
of quantity is rather involved. In the following we con- 
centrate on the contribution due to the ladder diagrams 
of the type shown in Fig. [6] which lead to an effective, 
screened interaction. These diagrams are connected at 
both ends as shown in Fig. [5] to obtain the noise. The 
lowest order contribution of this row is also contained 
in the perturbation theory as Fig. [2](4). The diagrams 
shown in Fig. [6] describe processes in which an electron 
pair is consecutively scattered on uncorrelated electron- 
hole pairs (bubbles). It is well-known that the neglected 
processes are unimportant in the case of bulk extended 
fermionic systems of high density and the respective ap- 
proximation is the familiar random phase approximation 
(RPA) . In the RPA approximation one sums over the di- 
agrams with leading contribution for each order. This ar- 
gument does not apply in our case as a zero-dimensional 
dot the concept of a density does not exist. However, 
with help of the GF in Eq. Q one can show that more 
convoluted diagrams should contribute less than the sim- 
ple ones we use in this approximation, at least for not too 
big ui. 

Furthermore, diagrams similar to (1) and (2) from 
Fig. [2] are composed of two single particle GF and do 
not spawn vertex corrections. The frequency-dependent 
noise however is an inherent vertex dependent quantity 
and effects from diagrams of this kind are expected to 
only result in features which are already well known. To 
support this point we calculated the noise also with help 
of a resummation using the analytical non-equilibrium 
self-energy^ which revealed no qualitative contribution. 
Moreover, we investigated diagrams related to (3) from 
Fig. [2] by using the screened interaction of Fig. [6] in dia- 
gram (3), Fig. [2} Numerical checks for various parameter 
constellations confirm that there is no significant contri- 
bution from these diagrams. This strongly suggests that 
the features we find using our resummation of diagrams 
similar to Fig. [2](4) are not canceled by diagrams which 
are neglected and are part of the exact solution of the 
model. 

Although our approximation is very similar in its spirit 
to the RPA, a precise identification of the small param- 
eter is impossible. Nonetheless, we can roughly estimate 
the relevant parameter regime by observing that in all di- 
agrams we take into account the loops describe electron- 
hole production/annihilation processes, in which every 
electron is annihilated by its 'own' hole. Since any re- 
combination process with a different ('alien') hole can 
only take place by tunneling into the electrode and back 
the recombination processes with 'own' holes are more 
probable for larger dwelling times r of the electrons on 
the dot. Since r ~ 1/T our approximation is valid when 
r is small in comparison to U . 

On the other hand, Fermi liquid theory for the non- 
equilibrium transport through an Anderson impurity * 11 ^ 23 ! 
implies that the principal observables can be expressed in 
terms of equilibrium susceptibilities. This is just as in the 
equilibrium case which was treated in the famous series of 



contributions by Yamada and YosidaP^H Interestingly, 
from the diagrammatic point of view the equilibrium sus- 
ceptibilities are represented by simple loops. Although it 
is very difficult to establish a one-to-one correspondence, 
our approach is a version of a Fermi liquid theory applied 
to non-equilibrium noise and it is therefore expected to 
yield adequate results particularly in the low energy sec- 
tor. 

The sum over ladder diagrams shown in Fig. [6]is equiv- 
alent to the solution of a set of 32 linear equations: 

r^ m "( Wl , W2 , W 3)« 

iUY J PDt P ^i)Dl\u 2 )T^ mn {u 1 ,u 2 ^), (12) 



r^™"^!,^,^) « 2n5(uj 3 -Lo 2 )D kn (u 1 )D™ l (u 3 ) 

+iU^2 pDo P (u>i)D p l (cj 2 )m mn ( Wl , W2, w 3 ) . (13) 
P =± 

An explicit analytical solution is possible but rather 
lengthy so we discuss the respective numerical results 
only. From general grounds, Y klmn should satisfy certain 
symmetry conditions. Similar to other partial resumma- 
tion schemes, it vanishes for uj 2 ^ u 3 . Moreover, it should 
satisfy T klmn (u} U u} 2 , uj 2 ) = T mnkl (uj 2 , lo u wi) due to the 
symmetry under exchange of the incident particles. Be- 
cause of time-reversal symmetry it also suffices the sym- 
metry when changing u>i <-> —Ldi, for i = 1,2. Last but 
not least, the zero- frequency noise in equilibrium has to 
vanish. 

The vertex function possesses simple poles when U > 
2T. This is due to the fact that we are using the cutoff- 
free wide flat-band model. The numerical evaluation of 
the noise according to the prescription of Eq. ( 10 ) there- 



fore needs regularization, which is achieved by imposing 
a hard cutoff scheme which models more realistic finite 
bandwidths. 
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FIG. 7: Interaction correction to the noise. The interaction 
strength is set to U = 0.99 where F ktmn is regular. The V = 
curve is enlarged by a factor of 10. 

We consider the two different regimes U <> 2T sepa- 
rately. In the U < 2T case where Y klmn is completely reg- 
ular, we obtain the results shown in Fig. [7]for the interac- 
tion induced correction to the non-interacting noise spec- 
trum. In equilibrium, finite interaction causes a small 
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suppression of the noise similar to the result of the sec- 
ond order perturbation theory. For larger voltages we 
observe an enhancement in the domain M < V. 
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FIG. 8: Interaction correction to the noise with U = 1.5 in 
the non-regular regime of Y klmn . The V — curve is enlarged 
by a factor of 10. 

More interesting is the case of strong interactions U > 
2T. As mentioned before, Y klmn develops simple poles 
in this regime that are also found in the integrand of 
Eq. (10) as long as ui < V/2. Due to these poles the 
correction to the non-interacting noise (Fig. [8]) develops 
sharp peaks at positions ±w ros < V/2 that do not seem to 
be affected by V. This turns out to be a generic feature 
and the distance between the peaks is solely defined by 
the interaction strength. The peaks disappear for too 
small voltages though, see Fig. [8j 
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FIG. 9: The full noise power for V = 4 and different interac- 
tion strengths calculated with the diagrammatic resummation 
scheme. For (7 = one recovers the step-structure. 

Figure [9] shows the full noise spectrum in the non- 
equilibrium configuration for different Coulomb interac- 
tion strengths. With growing U two distinct peaks at 
frequencies ±w res develop from the non-interacting step- 
like graph which are not seen in equilibrium (compare to 
Fig. [8]) . Again, the two peaks correspond to emission and 
absorption processes and the noise spectrum is known to 
exhibit similar peaks in other systems, see for example 
Ref. [3S] . The position of the peaks w res as a function of U 
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FIG. 10: Position of the peaks in the noise spectrum for 
V — 20 and various values of U. The black line shows a 
%/U behaviour. 



is highly non-analytic and follows a square-root law with 
high precision. This fact points towards the generation 
of a new energy scale (in a non-perturbative way because 
otherwise it would be algebraic in U), which appears to 
play a fundamental role for the dynamical properties of 
the system. The only plausible option is the Kondo tem- 
perature given by^ 



T K = exp (-7TU/8T) 



(14) 



is plotted in Fig. 10 Interestingly, the resulting behavior 



The presence of pronounced maxima in the non- 
equilibrium noise at frequencies ~ Tk can be well- 
understood. At this energy the system is very sensitive to 
fluctuations since the absorbed/emitted energy is com- 
parable to the formation energy of the Kondo singlet, 
thus making it an efficient energy dissipation channel. 
However, this effect is not directly related to other phe- 
nomena of the Anderson model as for example the zero 
bias conductance peak to which the results of Ref. [40] , 
which presents data in the different regime V <C U, are 
connected. 



IV. SUMMARY 

In conclusion, we have investigated the frequency- 
dependent noise in the non-equilibrium Anderson impu- 
rity model. Perturbative results in second order for small 
interaction strength U show that the noise is suppressed 
at positive frequencies in equilibrium, compared to the 
non-interacting setup. In the non-equilibrium situation, 
however, we observe an increase in the noise for frequen- 
cies related to the applied voltage. Using a diagram- 
matic resummation scheme we could gain knowledge in 
the parameter regime which is non-perturbative both in 
7 and U. This approximation shows again a suppression 
of noise in equilibrium but sharp peaks for finite V at 
positions that are related to the Kondo temperature in- 
stead of the voltage. This is a consequence of the Kondo 
temperature defining the energy scale for internal fluctu- 
ations of the system. 
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